arXiv: 1504.05203v2 [gr-qc] 1 May 2015 


Collapse of massive fields in anti-de Sitter spacetime 


Hirotada Okawa^’^ , Jorge C. Lopes^ , Vitor Cardoso^’^ 

^ Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto, 606-8502, Japan 
^ Advanced Research Institute for Science & Engineering, 

Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan 
® CENTRA, Departamento de Fisica, Instituto Superior Tecnico - 1ST, 
Universidade de Lisboa - UL, Avenida Rovisco Pais 1, 1049 Lisboa, PoHugal and 
Perimeter Institute for Theoretical Physics Waterloo, Ontario N2J 2W9, Canada 

Gravitational collapse in asymptotically anti-de Sitter spacetime has a rich but poorly-understood 
structure. There are strong indications that some families of initial data form “bound” states, which 
are regular everywhere, while other families seem to always collapse to black holes. Here, we in¬ 
vestigate the collapse of massive scalar fields in anti-de Sitter, with enlarged freedom in the initial 
data setup, such as several distinct wavepackets, gravitationally interacting with each other. Our 
results are fully consistent with previous findings in the literature; massive fields, which have a fully 
resonant spectra, collapse at (arbitrarily?) small amplitude for some classes of initial data, and form 
oscillating stars for others. We find evidence that initial data consisting on several wavepackets may 
allow efficient exchange of energy between them, and delay the collapse substantially, or avoid it al¬ 
together. When the AdS boundary is artificially changed so that the spectrum is no longer resonant, 
cascading to higher frequencies may still be present. Finally, we comment on the asymptotically 
flat counterparts. 
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I. INTRODUCTION 

Numerical evidence indicates that general classes of 
initial data in anti-de Sitter (AdS) spacetime collapse 
to black holes (BHs), at very (perhaps even arbitrar¬ 
ily) small amplitudes [I|. One crucial ingredient for this 
phenomena is the confinement property of the timelike 
boundary of global AdS. Although not rigorously proven, 
this is thought to be one necessary condition for collapse 
at arbitrarily small amplitudes. As some of us pointed 
out, the results in AdS geometries raise the intriguing 
possibility that analogous effects might play a role in our 
univers^ by triggering collapse of compact stars, for ex¬ 
ample [3. In fact, similar characteristics - collapse at 
very small amplitudes - were observed in a “box” in flat 
spacetime under different boundary conditions Bi- 

By using a perturbative expansion of the field equa¬ 
tions, the original study [l| suggested a second important 
ingredient: that a fully-resonant “ or commensurable - 
spectrum resulting in a secular term is a necessary con¬ 
dition for collapse at arbitrarily small amplitudes. Com- 
mensurability is a condition on the proper eigenfrequen- 
cies of the background spacetime. In particular, there ex¬ 
ist combinations of i,j, fc, / G Z to drive arbitrarily-small 
initial pulses away from the perturbative regime Bi 
such that the corresponding modes satisfy 

uji = ujk + - uJi, ( 1 ) 

This condition guarantees that higher frequency modes 
are excited at full nonlinear level, potentially triggering 
a turbulent cascade and collapse. Such cascade to higher 
frequencies was indeed observed in the original and subse¬ 
quent works Bl- Nevertheless, because commensurability 
is generically violated at nonlinear level, it would imply 


- should it be a necessary condition for collapse - that 
generic families of initial data do not collapse at small 
amplitudes. In line with this reasoning, there is indeed 
strong evidence that other classes of initial data do not 
collapse, forming instead bound states or stars BB- 

In an attempt to understand the importance of “turbu¬ 
lent” gravitational behavior in asymptotically flat space- 
time, some of us found a surprising behavior. In a toy 
model describing a scalar field confined in flat spacetime 
by a thin shell of matter, we found that collapse ensued at 
very small amplitudes and for a variety of boundary con¬ 
ditions. For generic boundary conditions, the commen¬ 
surability property is not satisfied, there were no hints 
of systematic cascading to higher frequencies, but we 
still found consistent approach to gravitational coll^se, 
something at odds with the current understanding (2|. 

The numerical implementation used in Ref. Bl differs 
from those of other works on the subject, and so does 
the specific setup. As such, our purpose here is to in¬ 
vestigate collapse of scalar helds in AdS, using similar 
numerical procedures to those reported in the literature, 
but allowing more freedom in the system so as to under¬ 
stand whether the conclusions in Ref. B "were particular 
to that setup. Our purpose is to explore the gravita¬ 
tional collapse of scalars in AdS, but allowing some more 
freedom in the theory and in the setup, in particular 

(i) the scalar is generically massive, 

(ii) the outer boundary may be located at a prescribed 
position, rather than only at spatial infinity (following 
Ref. B); and 

(iii) The initial data will consist, generically, on a super¬ 
position of gaussian-like wavepackets. 

Condition (i) still implies a fully resonant spectrum for 
the scalar fluctuations, while condition (ii) breaks it. 
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II. SETUP 


The system we are studying is a generalization of a 
system studied by other authors [I| . In our work we will 
use c = h = G = 1. The model considered in this study 
describes the dynamics of a spherically symmetric mas¬ 
sive real scalar field in a 3-1-1 dimensional AdS spacetime. 
In such case the equations that rule the dynamics of the 
system are the Einstein equation and the Klein-Gordon 
equation, with a mass term, 

Gab + Agab = StT {Va^pVbCp - - QabV{<!>)) , (2) 

-^=daV^9‘'^db(j) - ^ = 0. (3) 

We want to solve these equations for massive fields. As 
such we need to define what is their potential. We will 
set the massive potential to V{(j)) = 

In our study we will use as ansatz the line element 


;2 

ds^ = -^— (Ae~^^dt^ + A~^dx^ + sin^ x (4) 

COS"' X 


where = —3/A and is the differential solid angle of 
the unit two-sphere. We can recover the radial coordinate 
in Schwarzschild coordinates r = Itanx. The symmetric 
nature of the system demands that its variables. A, S and 
4>, depend only on the time and radial coordinates, t and 
x. 

Using the metric ansatz above together with the auxil¬ 
iary variables $ = </)' and 11 = A~^e^(l), the Klein-Gordon 
equation ((S]) yields 


$ = (Ae-^n)', 

n = —\— (tan^ xAe~^^) 
tan X 

The Einstein equations yield also 



cos^ X 


(5) 

( 6 ) 


A=2(l+.2^in^")(l_^) 

sin 2x 

— 47r sin x cos xA -|- 11^) — tan x, (7) 

d'= — 47rsina;cosx -I-n^) . (8) 


By appropriately re-defining /i one can re-scale the AdS 
radius I out of the problem. In other words, our results 
scale completely with 1. When /x —>■ 0 one recovers previ¬ 
ous system of equations [l|. 

Eor the solution to be physically meaningful, the pro¬ 
files must be smooth everywhere. At the origin we set 
the boundary conditions 


(j){t,x) = +0(x^), 

(9) 

A{t, x) = 1 -1- 0{x^), 

(10) 

S(t, x) = 0 -1- 0{x^) , 

(11) 

and, at spatial infinity, we require that 

~ (/)oo(t)/ -bO(/+^), 

(12) 

A{t,p)c^l-Mp^ + Oip^), 

(13) 

s{t,p) ~ (5oo(t) + e>(p^'=), 

(14) 


where p=^—x and + The behavior 

for the scalar field is dictated by regularity of the Klein- 
Gordon equation in an asymptotically AdS spacetime. 


A. The eigenfrequencies 


To compute the eigenfrequencies of scalar fields in AdS 
at linearized level, it proves easier to work with AdS in 
(re-scaled, with the AdS radius I = 1) global coordinates, 

ds^ = —(r^ -I- l)dt'^ -\— ^ -h r'^dfl'^ . (15) 

r"' -I- 1 

The spherically symmetric, linearized Klein-Gordon 
equation for a field (j> = jr reads 


/ (2 + M^)) - (16) 

with / = r^ + l. The above equation is of hypergeometric 
type. The solution which vanishes at the origin is 


v]/ = r(r^ -1-1) 2 i? 


3 1 

1 

CO 

k 

UJ 

3 

to 

2 

2’ 
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, (17) 


Imposing Dirichlet conditions at infinity one finds 


ujl = k + 2n , 


(18) 


so that the spectrum is still fully resonant in the ter¬ 
minology introduced in the Introduction, and will give 
rise to secular terms 0. One can break the resonant 
spectrum by imposing Dirichlet conditions instead at a 
finite radius, as done in Ref. Q- Examples of resonant 
frequencies are shown in Table HI 


B. Initial data 

In this paper, we use two types of initial data. One 
is the single wavepacket initial data used in the litera¬ 
ture [H, 

$(0, s) = 0, n(0, x) = —e"t‘C 2 " , (19) 

TT 


TABLE I. Resonant frequecies for massive scalars in AdS, 
when the boundary is truncted at a finite rmax/Z = tanXmax. 


3^niax/^ 

UJ 1 


pi = 0 

pi = 2 


3.0000 

4.0000 

7r/2 

5.0000 

6.0000 


7.0000 

8.0000 


3.6497 

4.4463 

0.75 

7.1619 

7.7302 


11.056 

11.471 


9.3352 

10.198 

0.25 

20.452 

21.145 


32.455 

32.964 
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where a represents the width and e the wavepacket’s am¬ 
plitude. The other is composed of several wavepackets, 


^ 9 

n(0,a:) = ey^ exp 

i=l 

where we still assume $(0,x) = 0 and, the width and 
location of the i — th wavepacket are defined by ai,ri. 
We also define overall amplitude by e and the ratio of 
amplitude to the first wavepacket by Ci, with ai = 1. 


—4(tana; — 


( 20 ) 


III. NUMERICAL RESULTS 
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A description of the methods used to solve the astern 
in the massless case are available in the literature . We 
generalized the approach, which uses a RK4 time evolu¬ 
tion method with finite differences to compute the spatial 
profiles at each step of the evolution. Note that it is pos- 
sile to identify the apparent horizon (AH) formation at 
XAH corresponding to a zero of A(t, x) or a point at which 
A{t,XAH) falls below a prescribed thresholdfl^. Gener- 
ically, our results show second-order convergence for the 
evolution of massive fields (and fourth-order for massless 
fields). A convergence test is shown in Appendix 


A. Massive field collapse 
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Our results for the gravitational collapse of a sin¬ 
gle wavepacket of massive fields are best summarized 
by Fig. [TJ Here, we compare directly the BH radius 
formed through the collapse of a massless(/rZ = 0) and 
a massive(/r/ = 2) scalar for a relatively small width 
a = //16 (top plot) and large amplitude e/ = 4 (bot¬ 
tom plot). For these, and despite the relatively large 
value for pi - larger than any other scale in the prob¬ 
lem - collapse proceeds almost blindly with respect to 
the field mass. Although not shown here, other mass 
terms (we tested for example pil = 1) also exhibit the 
same trend. The time elapsed until collapse also shows 
a remarkable agreement for massive and massless fields. 
There are some fine prints, visible in the lower panels of 
the plots. For fixed width and large initial amplitudes 
(within the “Choptuik” stage of prompt collapse), the 
mass term tends to make the BH size larger. However, 
the BH seems to become smaller almost everywhere else 
at smaller amplitudes, with the exception of the critical 
points. 

B. Stability islands 

Stability “islands” corresponding to bound states, or 
boson stars in AdS, where no gravitational collapse is 
observed, were previously reported for massless fields @ . 
To investigate these possible end-states we have varied 
the width of the initial scalar profile, and monitored the 


FIG. 1. Direct comparison by BH radii formed through 
massless collapse(/il = 0) and massive collapse(/rl = 2). 
The width/amplitude of the Gaussian wavepacket is fixed at 
cr = 1/16 (top) and el = 4 (bottom). 

collapse time. The results are summarized in Figs. [5]13] 
for massless (/il = 0, upper panels) and massive {pi = 2, 
lower panels) fields. The results for massless and mas¬ 
sive fields differ quantitatively, but not qualitatively. In 
particular, we still find strong indications that there is 
no collapse for families of initial data with moderately 
large width. For example, a = 21/5 takes an extremely 
long time to collapse below an amplitude ecr ~ 0.65; on 
the other hand for a massive field with pi = 2, not only 
a = 21/5 but also cr = 1/4 show extremely delayed col¬ 
lapse below ecr ~ 0.62. There is thus a transition to 
new nonlinear solutions of the field equations; these are 
known as oscillatons in asymptotically flat space m and 
have been observed to form in the collapse of massive 
real fields [H. In asymptotically AdS geometries oscil¬ 
latons or “scalar field breathers” exist even for massless 
fields [l^ and our results show that they are generically 
approached for some families of initial data. The ap¬ 
proach of collapsing configurations to oscillatons is not 
completely trivial as some of these oscillate for a long 
time before collapsing, presumably because they are ei¬ 
ther linearly or nonlinearly unstable [T^ . This trans¬ 
lates into a nontrivial collapse time, as the amplitude 
decreases, visible in the inset of Fig. [21 
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FIG. 2. Critical behavior in AdS with different initial pulse 
width (T = //16, Z/8, Z/4 and 2//5. Upper panels denote the 
BH radius while lower panels show the collapse time as a 
function of initial amplitude. Top: massless fields, fj,l = 0. 
Bottom: massive fields, fj.1 — 2. 

It should also be noted that one can find whether such 
bound-states arise or not by decreasing the initial ampli¬ 
tude and monitoring the final BH size near the critical 
amplitude, as in Fig. [21 Our results imply that the tran¬ 
sition from BHs to boson stars or bound-states is similar 
to that reported for the collapse of massive fields in a 
Minkowski background [HI, [T^ : the BH near the criti¬ 
cal amplitude has a finite mass. The confining box size 
is now the AdS boundary, while in previous Minkowski- 
background studies, it was simply the wavepacket width. 

Moderately large-width initial data does not seem to 
collapse, and bound-states arise. As pointed out in 
Ref. [i^, much larger widths actually restore the col¬ 
lapse, this is apparent in Fig. [31 for massless fields. Cor¬ 
responding initial data are also shown in the lower panel. 
Although we do not show a detailed mode analysis, the 
transfer of energy to higher-frequency modes for “nor¬ 
mal” (i.e., as reported originally in Ref. Q) = V4) 
and “restored” {a = Al) collapse, is already clear from 
the evolution of the Ricci scalar time evolution at the 
origin. By contrast, the “delayed collapse” configuration 
shows nontrivial oscillation patterns. We note that the 
initial data for restored collapse can be regarded as a 



t/i 



x/l 

FIG. 3. Top Panels: Ricci scalar evolution at the origin for 
“normal” (top), “delayed” (middle) and “restored” collapse 
(bottom) of massless fields. For very large widths, we find 
that collapse to BHs ensues again and that energy cascade to 
higher frequencies occurs efficiently. Bottom Panels: Gorre- 
sponding initial profiles of the scalar field H. 

localized wavepacket with respect to the corresponding 
density distribution shown in Ref. [l^ . 

For a region where delayed collapse is found, we fitted 
the critical amplitude Cc which divides bound-states from 
BHs in Fig. 01 We find that Cc ~ Finally, the 

bottom panel shows that in a re-scaled axis plot, restored 
collapse (cr = 41) exhibits a completely different behavior 
from “normal” (cr = I/A) or delayed collapse(cr = 1). 

C. Finite boundary and resonant frequencies 

The behavior reported previously, where some classes 
of (small-width) initial data collapse at very small am¬ 
plitudes whereas others (larger-width) do not, is not de¬ 
pendent on the resonant spectrum of AdS. This can be 
tested by introducing an artificial “mirror” in the AdS 
bulk, which amounts to imposing Neumann boundary 
conditions for (f) at that location. The collapse properties 
are summarized in Figs.jSJTl 

In Fig. O we show time evolutions of the Ricci scalar 
at the origin for a massless(^/ = 0) and massive)^/ = 2) 
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FIG. 4. Delayed collapse of massless fields: Critical am¬ 
plitude for collapse as a function of width (top). The bot¬ 
tom panels show the collapse time and BH size for re-scaled 
widths, where the contrast between small (cr = //4), interme¬ 
diate {a = 1) and very-large {a = 4/) initial width is apparent. 
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FIG. 5. Time evolution of the Ricci scalar 11^ at the origin 
during massless(red solid line) and massive(blue dotted line) 
collapse with an artificial boundary at Xma^/l = 0.75. Am¬ 
plitude and width of the initial pulse are set by el = 6.5 and 
O’ = //16, respectively. 
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collapse with an artificial mirror boundary at Xmax/l = 
0.75. Notice how the scalar curvature increase super- 
exponentially for both cases, despite the fact that both 
clearly have a non-resonant spectra. 

The structure of collapse time as function of ampli¬ 
tude and width shows an interesting structure depicted 
in Figs. inilTl For a wide range of initial conditions, the 
time development is not strongly dependent on a mass 
term, as is evident from Fig. [SI We also show the rescaled 
collapse time with an artificial mirror at Xmax in Fig- [3 
Delayed collapse is clearly seen for Xmaxll = 0.25, while 
the collapse time with finite boundary for a relatively 
large Xmax also scales as e“^, for example, Xmaxll = 0.75. 

In summary, there are certainly families of initial data 
which do not collapse, in agreement with the view that 
backgrounds with non-resonant frequencies cannot cause 
collapse at arbitrarily small frequencies. On the other 
hand, we do find collapse for other generic families of ini¬ 
tial data at finite amplitude (as must be the case numer¬ 
ically). In addition, we continue to observe migration to 
higher frequencies and a collapse time which is of the or¬ 
der of the full (resonant) AdS spacetime. Unfortunately, 
our results do not shed new light on this issue, and the 


FIG. 6. Critical collapse with a finite boundary for mass- 
less(top, pi = 0) and massive (bottom, pZ = 2) scalars. The 
outer boundary is located at a prescribed point, by intro¬ 
ducing an “artificial mirror” at Xmax- The width of scalar 
wavepackets is set hy a = 1/16. For Xmax/l ~ 7r/2 one recov¬ 
ers Fig. (2] 

results in Refs. continue to be poorly understood 

(unless one admits that collapse will eventually occur at 
a very small but finite amplitude, whose threshold is at 
the moment unpredicted). 


D. Energy transfer between wavepackets 

Our results do shed some more light on the existence 
of nonlinearly stable bound states, and the role of grav¬ 
itational redshift in their formation. For this purpose, 
we now focus on the collapse of initial data described by 
more than one wavepacket. 

In Fig. [51 for example, we show the time it takes for col¬ 
lapse, starting from initial data described by two (upper 
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FIG. 7. Rescaled collapse time for massless fields with a width 
a = Z/16. Different marks correspond to different mirror loca¬ 
tions as Fig. |6l Dotted line was obtained by fitting the results 
from Xmaxjl ~ 0.75 simulations. 
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FIG. 8. Collapse time for two different initial data sets: two 
(top) and three-wavepacket initial data (bottom). The col¬ 
lapse time is shown as a function of the overall amplitude el 
in (Uni). Parameters are set by = {Z/8 ,/tt/S, 1} 

and {<J 2 ,r 2 ,a 2 } = for two wavepackets and by 

{o-i,ri,ai} = {Z/4, 0,1}, {cr 2 ,r 2 ,a 2 } = {//4, / tan(7r/8), 1/4} 
and {cr 3 ,r 3 ,a 3 } = {Z/4, Z tan(7r/4), 1/8} for three wavepack¬ 
ets, respectively. Crosses (Circles) denote the results with 
low(high) resolution. Notice that, although the three- 
wavepacket initial data does not collapse, the width of each 
individual wavepacket is very small, and would lead to col¬ 
lapse were they in isolation. 

panel) and three wavepackets (lower panel). It is appar¬ 
ent from these plots that even relatively large amplitudes 
need some more reflections to collapse, as compared with 
a single wavepacket. In fact, the three-wavepacket config¬ 
uration clearly displays delayed collapse, and exhibits a 
threshold amplitude below which no black hole is formed. 
We should also add that the importance and effect of 
wavepacket superposition for gravitational collapse are 
also discussed in Ref. (l^ . 

Typical evolutions of Ricci scalar monitored at the ori¬ 
gin are exemplified in the left column of Fig. [5] and corre¬ 


FIG. 9. Ricci scalar evolution and corresponding initial data. 
Parameters are the same as in Fig. |8l and the overall ampli¬ 
tudes are given by el = 1/5 for two pulses and el = 1 for 
three pulses. In the top panel, collapse begins in the smaller- 
amplitude wavepacket. 



0 71/8 71/4 371/8 71/2 


x/l 

FIG. 10. Early phase evolution of Ricci scalar at the ori- 
gin(top) using three wavepackets initial data and snapshot of 
lapse function at t/l = 0.0,8.4,16.7,24.5 (bottom). Parame¬ 
ters are the same as those in Fig. 


sponding initial data are shown in the right column. The 
Ricci scalar grows monotonically for two wavepackets, 
while it oscillates almost periodically for three wavepack¬ 
ets. A mode analysis of the three-wavepacket evolution 
shows a cascade to higher frequencies at early times, and 
a constant frequency content at late times. 

In order to understand the finer structure of the evolu¬ 
tion, the top panel of Fig. [TU] shows only the early phase 
evolution of Ricci scalar. There are some features worth 
highlighting. 

(i) The initial (three-wavepacket) configuration corre¬ 
sponds roughly to the three peaks observed before t/l <\ 
and is similar to that within < t/l < A (after one re¬ 
flection) and 6 < t/l < 7 (after two reflections). How¬ 
ever, the initial configuration produced, besides these in¬ 
going waves, out-going waves. These are observed around 
2 < t/l < 3 or at around t/l = 5 and correspond to 
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waves directed to the AdS boundary. Let us call them 
“secondary waves”. 

(ii) The original in-going packets (say at t = 
0,3, 6, 8.5,11.5) are always accompanied by high am¬ 
plitude peaks to their right. The bottom panel of Fig. [TUI 
shows that these are in a higher gravitational poten¬ 
tial, and can work as a source of redshift for those 
packets. Conversely, these high-amplitude peaks can be 
blueshifted. In other words, these localized packets can 
become broader or sharper, depending on the gravita¬ 
tional potential that they encounter. 

(iii) However, we find that the secondary waves can “mi¬ 
grate” and even take-over the originally in-going pulses 
(in this example, after t ^ 17). As a result, the red- 
shift/blueshift pattern is inverted and processes contin¬ 
ues in a cyclic way. 

Since the scalar field is responsible for the spacetime 
curvature in this setup, general-relativistic effects should 
become important when high-amplitude waves are lo¬ 
cated near the origin. The bottom panel of Fig. [TUI 
shows the lapse function defined by a = \fAe~^ j cosx 
in the line element ®. What this really shows is how 
clocks tick differently as a function of coordinate x for 
different instants, and can explain the migration of sec¬ 
ondary waves as mentioned above, considering that high- 
amplitude waves act as a source of gravitational poten¬ 
tial and other waves apart farther from the origin move 
faster than the ones near the origin. Also, the time 
dependency of lapse profile gives evolutions of original 
and secondary waves, because the frequency of ingo¬ 
ing (outgoing) waves in the gravitational potential can 
be changed by gravitational blueshift (redshift). For a 
free-fall observer, the difference of frequencies between 
waves dX X = Xi and a: = 0:2 is described formally by 
W 2 /W 1 = a{xi)/a{x 2 ) if the background spacetime is 
fixed. That implies that higher-amplitude and wider- 
width waves could cause the frequency shift much more 
efficiently than the ones slightly farther from the origin. 
The existence of many such wavepackets would give a 
meta-stable state, continuously exchanging energy. 

The above is, of course, a tempting explanation of why 
this type of initial data does not lead to collapse, but is 
far from rigorous. 


IV. DISCUSSION 

We investigated the collapse of spherically symmetric 
scalar fields in AdS for massive fields for generic initial 
data, and including an artificial “mirror,” extending pre¬ 
vious studies on the subject. Our main motivation was 
to understand role of “commensurability” in the process 
and the mechanism leading to collapse or to stable con¬ 
figurations. Our main findings are that 

I. In parallel with previous studies, we also find that 
there are stability islands for some regions of parameters 
in the space of initial conditions. In these regions - both 
for massless and massive fields - we observe quasi-stable 


bound states, or oscillating configurations. 

II. Single-wavepacket initial data for moderately-large 
widths evolve into nonlinearly-stable configurations at 
very small amplitudes (see also Fig. 7 in Ref. @) and 
at finite amplitude it results in a delayed collapse, as 
summarized in Fig. [T] However, we also confirm that 
much wider initial data restores the collapse again, in¬ 
dicating that there are thresholds dividing the phase of 
nonlinear solutions into two final states. Unlike turbulent 
collapse, the threhsold amplitude yields a finite-mass BH, 
in a behavior very similar to that found in the collapse of 
massive fields in asymptotically flat spacetimes [IJ, [ij . 

HI. Even if the condition of commensurability is vi¬ 
olated, for example by inserting an artificial mirror at 
a finite distance, our results exhibit energy transfer to 
higher modes and indicate the “usual” t ^ depen¬ 
dence (see Fig. El). In addition, for smaller mirror radius, 
we again observe nonlinearly stable, oscillating bound 
states. Our results suggest that the threshold of the ini¬ 
tial pulse width for transitions is around (jjxraax ~ 0.2 
for single wavepacket initial data. 

IV. We have argued that other, more generic type 
of initial data can produce stable bound states by 
an interplay of blueshift/redshift effects, and that this 
mechanism may even be at play in large-width, single- 
wavepacket data. For example, data consisting on two 
wavepackets seem to efficiently transfer energy from one 
wavepacket to the other, significantly delaying their col¬ 
lapse as compared to the single-wavepacket data. Inter¬ 
estingly, three-wavepacket initial data can evolve into sta¬ 
ble bound-states, very much like moderately large-width 
initial data consisting on a single wavepacket. 

We should highlight the fact that the collapse of mas¬ 
sive field exhibits qualitatively the same behavior as that 
of massless field in our mass parameter choice; it is pos¬ 
sible that much larger field mass produces more easily 
bound states, which would be the analog of flat-space os- 
cillatons. As a final comment, we should stress that our 
findings, both in this work and in previous Q, flat-space 
studies are perfectly consistent with the prediction that 
gravitational collapse at arbitrarily small amplitudes re¬ 
quires commensurable spectra Non-commensurable 
spectra yields collapse at finite amplitudes, which might 
be of some physical consequence in a variety of scenar¬ 
ios. An important open question concerns the threshold 
amplitude when the spectra is non-commensurable. 


Appendix A: Convergence test 

In Fig.[TTl we show convergence tests of the scalar field 
^ as a function of time for massive collapse corresponding 
to Fig. El Our results are compatible with second order 
convergence for (j) (obtained by a trapezoidal integration 
of <I>, cf. main text above ©)• 
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FIG. 11. Convergence tests for the scalar field (j) setting the 
parameters of initial pulse by ct = 1/16, = 2 and el = 6.8. 

The second order convergence with different resolutions Ni = 
N/2, N,2N is shown by defining the L2-norm, |0jvi,jv2(t = 

G)| = (0iv,(GO) - 
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